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Abstract 

The  transient  behavior  of  a  proton  exchange  membrane  fuel  cell  (PEMFC)  with  porosity  is  investigated  in  this  study  using  a  two-phase,  half-cell 
model.  The  thin  film  agglomerate  approach  is  used  to  model  the  catalyst  layer.  Both  vapor  transport  and  liquid  water  transport  in  the  PEMFC  are 
examined  in  this  study.  Proton  transport  is  much  faster  than  the  gaseous  and  liquid  water  transport.  The  ionic  potential  reaches  a  steady  state  level  in 
~10_1  s  but  liquid  water  transport  takes  ~10  s.  The  variation  of  the  ionic  potential  loss  reaches  a  critical  value,  decreasing  to  a  steady  state,  and  is  not 
monotonic.  The  gas  diffusion  layer  (GDL)  and  the  catalyst  layer  (CL)  porosity,  which  can  affect  cell  performance,  have  been  carefully  investigated. 
The  current  density  rises  rapidly  within  10-2  s,  then  remaining  constant.  After  1  s,  this  is  affected  by  the  cell  voltage,  GDL  porosity,  and  CL  porosity, 
and  if  the  GDL  porosity  is  below  0.4,  the  current  density  drops.  For  the  gas  diffusion  layer  porosity,  the  current  density  increases  between  £GDL  =  0.2 
and  £qdl  =  0.5,  with  increased  GDL  porosity.  For  the  catalyst  layer  porosity,  the  optimum  value  appears  between  £Cl  =  0.06  and  £Gl  =  0.1. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

In  recent  years,  fuel  cells  have  been  actively  developed 
for  use  in  portable  devices  including  mobile  phones,  com¬ 
puter  notebooks,  power  tools,  digital  cameras,  etc.,  proton 
exchange  membrane  (PEM)  fuel  cells  have  many  advantages, 
such  as  a  potentially  lower  cost  per  kW,  fast  start-up,  and  a 
lower  operating  temperature.  These  applications  require  a  highly 
dynamic  model,  but  most  recent  fuel  cell  studies  use  steady  state, 
rather  than  transient  models.  Developing  a  transient  model,  for 
obtaining  higher  PEM  fuel  cell  performance  is  therefore  very 
important. 

A  PEMFC  is  a  sandwich-like  structure  consisting  of  a  mem¬ 
brane,  gas  diffusion  layers  (GDL),  and  catalyst  layers  (CL).  The 
catalyst  layer  is  more  complex  than  the  other  layers,  because  it 
is  here  that  electrochemical  reactions  take  place,  and  different 
phase  types  are  present.  There  are  many  approaches  for  study¬ 
ing  the  catalyst  layer.  If  it  is  considered  to  be  an  interface  [1-3], 
many  parameters  will  be  ignored,  and  the  results  will  be  inac¬ 
curate.  Generally  speaking,  the  results  will  be  higher  than  they 
should  be.  Various  approaches  have  been  proposed  to  counter 
this,  such  as  the  thin  film  model  [4]  and  agglomerate  model 
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[5-10].  In  the  thin-film  model,  catalyst  particles  are  covered  by 
a  polymer  electrolyte  film,  and  gas  pores  exist  only  within  the 
electrode.  Film  thickness  is  uniform  and  very  small  compared 
to  pore  size.  In  the  agglomerate  model,  catalyst  particles,  elec¬ 
trolytes,  and  gas  pores  form  a  homogeneous  mixture.  Several 
researchers  [11-14]  compared  their  model  to  experimental  data. 
Broka  and  Ekdunge  [11]  have  suggested  that  the  agglomerate 
model  is  more  accurate  relative  to  other  models. 

Most  models  treat  the  catalyst  layer  as  a  single  phase.  It  is 
in  fact,  multiphase  in  a  fuel  cell.  More  sophisticated  treatments 
have  been  used  for  a  two-phase  flow  in  fuel  cells  [15-19]. 
Van  Nguyen  and  co-workers  [16-19]  assumed  that  liquid 
water  transport  through  the  porous  electrode  is  driven  by  gas 
flow  shear  and  capillary  forces,  and  they  also  assumed  a  net 
flux  through  the  membrane.  Permeability  was  taken  to  be  a 
linear  saturation  level  function.  According  to  the  two-phase 
model,  if  liquid  water  is  neglected,  the  cell  performance  will 
be  overestimated.  Portable  device  applications  require  highly 
dynamic  loadings  so  PEMFC  steady  state  models  must  have 
some  imperfections.  Most  studies  take  a  single  factor  into 
consideration,  leading  to  overestimations. 

A  transient,  two-phase,  one-dimensional,  isothermal  and 
isobaric  model  of  a  proton  exchange  membrane  fuel  cell 
cathode  is  the  basis  for  this  study.  The  approach  offered  by 
Springer  et  al.  [3,20]  is  assumed  for  the  membrane.  The  thin 
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Nomenclature 

a 

surface  area  per  unit  volume 

b 

tafel  slope 

Cj 

concentration  of  species  j  (mol  cm-3) 

Dj 

diffusion  coefficient  of  species  j  (cm2  s-1) 

F 

Faraday  constant 

H 

Henry  constant  (atm  cm3  mol-1) 

i 

current  density  (A  cm-2) 

kj 

oxygen  reaction  rate  constant  (s-1) 

kc 

condensation  rate 

ky 

evaporation  rate 

K 

permeability  (cm2) 

M 

molecular  weight 

Nj 

molar  flux  of  species  j  (mol  cm-2  s-1) 

P 

pressure  (atm) 

r 

radius  (cm) 

s 

liquid  water  saturation  level  in  porous  medium 

T 

temperature  (K) 

Vs 

cell  potential  (V) 

y 

mole  fraction 

Greek 

thickness  (cm) 

£ 

porosity  or  volumetric  fraction 

0 

ionic  potential  (V) 

<P 

Thiele  modulus 

K 

conductivity  (£2-1  cm-1) 

X 

water  content 

F 

viscosity  (gem-1  s-1) 

P 

density  (g  cm-3) 

X 

tortuosity 

$ 

effectiveness  factor 

Subscripts  and  superscripts 

agg 

agglomerate 

CL 

catalyst  layer 

eff 

effective 

g 

gas  phase 

GDL 

gas  diffusion  layer 

MEM 

membrane 

N 

Nation  phase 

02 

oxygen 

p 

proton 

Pt 

platinum 

sat 

saturation 

V 

vapor  water 

w 

liquid  water 

film-agglomerate  approach  presented  by  Lin  et  al.  [19]  is 
applied  to  describe  the  catalyst  layer.  In  terms  of  modeling 
two-phase  transport  in  both  the  gas  diffusion  layer  and  catalyst 
layer,  the  approach  developed  by  Natarajan  and  Van  Nguyen 
[18]  is  applied.  Presenting  interactions  between  each  parameter 
and  the  cell  performance  in  the  transient  state,  such  as  gas 
diffusion  layer  porosity  and  catalyst  layer  porosity,  is  the 


objective  of  this  study.  It  is  hoped  that  appropriate  parameters 
will  help  generate  a  higher  performance  fuel  cell. 

2.  Mathematical  modeling 

This  study  is  based  on  the  model  of  Lin  et  al.  [19]  to  develop  a 
half-cell  transient  model  (Fig.  1)  and  investigates  how  the  GDL 
and  CL  porosities  affect  the  fuel  cell  performance.  It  also  ana¬ 
lyzes  electrochemical  kinetics  and  transport  of  oxygen,  vapor 
water,  and  liquid  water  in  the  GDL,  the  catalyst  layer,  and  in 
the  proton  exchange  membrane.  After  PEMFC  start-up,  air  dif¬ 
fuses  from  the  channel  to  the  catalyst  layer  via  the  gas  diffusion 
layer.  In  the  meantime,  oxygen  dissolves  into  the  Nation  film 
and  reaches  the  pellet  surface;  then  liquid  water  is  generated. 
The  following  assumptions  are  made: 

(1)  The  gas  phase  obeys  the  ideal  gas  law. 

(2)  In  the  GDL  and  CL,  electronic  resistance  is  negligible. 

(3)  Catalyst  pellets  are  treated  as  cylindrical  and  consist  of 
carbon- supported  platinum  and  Nation,  covered  by  a  Nation 
film. 

(4)  Catalyst  pellet  radius  and  Nation  film  thickness  are  uniform 
in  the  catalyst  layer. 

(5)  Catalyst  pellets  are  treated  as  homogeneous  and  oxygen  dif¬ 
fuses  into  pellets  via  Nation  film. 

(6)  The  Thiele  modulus  approach  is  applied. 

(7)  The  oxygen  reduction  reaction  generates  liquid  phase  water. 

(8)  The  ionic  potential  between  the  anode  catalyst  layer  and 
the  membrane  is  approximately  zero  because  of  the  fast 
hydrogen  oxidation  reaction  rate. 

Conservation  equations  are  listed  in  Table  1  with  five 
variables. 

(1)  The  concentration  of  oxygen  in  the  gas  phase  [Cq  ]. 

(2)  The  concentration  of  vapor  water  in  the  gas  phase  [Cf  ] . 

(3)  The  liquid  water  saturation  level  [ s ]  (the  ratio  between 
liquid  water  volume  and  total  void  volume  in  the  porous 
medium). 

(4)  The  concentration  of  liquid  water  in  Nation  phase  [C^J] . 

(5)  The  ionic  potential  [0]  (potential  in  Nation  phase). 

The  flux  expression  of  gas  phase  is 
Nj  =  -Dj[si0(l-s)]rVCgj  (1) 

£q  is  intrinsic  porosity  of  GDL  or  CL.  With  diffusion  coefficient, 
Dj ,  nothing  but  temperature  dependence  is  considered  [21].  The 
Bruggeman  correlation  is  used  to  assess  porosity  and  tortuosity 
effects. 

Water  interfacial  transfer  rate  between  liquid  and  vapor,  Rw , 
is  the  same  as  that  used  by  Lin  et  al.  [19].  Ro2  is  the  reaction 
rate  of  oxygen  in  the  catalyst  layer. 

-  = _ ElUoi _ c§ 

2  (8N/arD^2)  +  (Sw/arD^2)(H^2/Hg2)  +  (l/^kT)  °2 

(2) 

Hq2  and  Hq2  are  the  Henry’s  constants.  Hqo  is  the  oxygen 
between  air  and  the  Nation  phase  and  Hq2  is  the  oxygen  between 
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Fig.  1.  Schematic  of  the  model  domain  [19]. 


air  and  liquid  water.  8 n  is  the  thickness  of  the  Nation  film.  Dqo 
is  the  diffusivity  of  oxygen  in  the  Nation,  and  Dqo  is  the  diffu- 
sivity  of  oxygen  in  liquid  water.  ar  is  the  outer  surface  area  of 
agglomerates  per  catalyst  layer  unit  volume. 


ar 


2 

r Agg  +  <$n 


(l 


rAgg  is  the  radius  of  a  catalyst  pellet.  Liquid  water  forms  a  film 
around  the  top  of  Nation  film.  If  it  is  not  evaporated,  its  thickness 
can  be  estimated  by 


ar 


The  reaction  rate  constant 


pCL 
£0  * 


1 


4FC0 


agg  • 

a?bbio  exp 


2,eff 


2.3Q3(ys  -  0  -  Uref) 
b 

(5) 


The  active  catalyst  surface  area  per  unit  volume  of  agglom¬ 
erates 


agg  Opt^Pt 

Pt  =  5cl(1  -  ^L) 


apt  is  the  surface  area  per  mass.  mpt  is  the  catalyst  loading  of 
the  electrode.  <$cl  *s  the  thickness  of  catalyst  layer  and  b  is  the 
Tafel  slope. 


Table  1 

Governing  equation 


Variables  GDL 


CL 


MEM 


s 


0 


|((1  -s)«oDLco2)  = 

£>o2s^DLr[[(l  -  ^)]tV2C®2  +  V[(l  -  ^)]tVC«2] 
|((1-S)S°DLC®)  = 

-Rw  +  OvfigDLt[[(l  -  s)]rV2C®  +  V[(l  -  s)]rVC®] 

t  I  =  &{-*-£)  ^2s + (Vs>2>  +  «w 


gdl  _  nN\72rN 

m  a f  —  v  L-w 


a,((l  s)£o  Cq2)  —  R0l  + 

Do2sgr[[(l  -^)]rV2C*2  +  V[(l  -s)]rVC®2] 


|((l-S)£gcf)  = 

-Rw  +  Dv^Lr[[(l  -  S)]TV2C®  +  V[(l  -  5)]rVC®] 


pCL 

e0  ds 


M„  f  -  ,4"d'  +  2)«o2  +  Rw  + 

^f(-^)(SV2s  +  (V,)2) 


XL  9CW 


0  —  /)NV2CN 

in  dt  w  v 

KN,effV20  -  4FRq2  =  0 


CS2=0 
cf  =0 


5  =  0 


r  _  r)Ny72/HS  i  MdArNv72T/ 

£m  ~dT  ~  V  Cw  d - 7T“  V  K+ 

v20  =  o 


Table  2 


Boundary  conditions 


Variables 

v=o 

GDL/CL 

CL/MEM 

X  =  L 

Cs 

co2 

U02  “  uo2 

Vgo2 

CT 

,  Igdl  =  A^0o  |gdl 

<ICL=0 

cs  — 

0 

Cs 

f* g  /^-air 

V-V  Vy 

N*\ 

g 

gdl  =  Nw  |gdl 

WylCL  =  0 

cs  — 

co2 

0 

s 

s  =  0 

Nw\ 

GDL  -  A^wIgdL 

A^wIcl  —  0 

s  =  0 

r'  N 

CN  =  0 

KJ 

VN 

2VW 

Igdl  —  A0J|gdl 

<|CL  =  <|CL 

G'N  I  _ 

e  w  |mem  - 

_  XV, eq 
w(anodc) 

0 

— 

fplGDL  —  0 

ZpICL  =  ip  I  CL 

0  =  0 
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The  effectiveness  factor 
1  3  <p  coth  (3  (p)  —  1 
<p  3cp 

The  Thiele  modulus 


rAgg 


kr/O  ~  £nL) 


0 


D 


N 

02,eff 


In  the  anode,  protons  are  produced  and  transported  through 
the  membrane  to  the  catalyst  layer.  The  flux  expression  of 
charges  in  the  membrane  and  catalyst  layer  is 


*'p  —  K"N,eff^0 


Liquid  water  flow  in  the  porous  media  obeys  Darcy’s  law.  Its 
driving  force  is  only  the  capillary  force  and  the  flow  is  similar 
to  that  used  by  Natarajan  and  Van  Nguyen  [18].  Liquid  water 
permeability  dependence  on  the  saturation  level  is  linear. 


Pw  Kw,0 

MwfJw 


(10) 


pw,  Mw  and  /zw  are  the  density,  molecular  weight,  and  viscosity 
of  liquid  water,  respectively.  Kw,o  is  the  permeability  of  liquid 
water  at  100%  saturation  level.  —  (dpc/ds)  is  treated  as  a  constant, 
but  the  values  used  at  GDL  and  CL  are  different. 

In  the  GDL  and  CL,  liquid  water  diffuses  in  the  Nation  phase. 
In  the  membrane,  liquid  water  transport  is  due  to  the  water  con¬ 
centration  gradient  and  electro-osmotic  drag  and  similar  to  that 
considered  by  Lin  et  al.  [19]: 


Nn  = 

1  v  w 


n 


CL 


F 


i  _  nNvrN 

lp  v  ^ w 


(11) 


Initially,  all  variables  are  zero  except  for  the  liquid  water 
concentration  in  the  Nation  phase,  which  is  assumed  to  be  in 
equilibrium  with  vapor  water  in  the  gas  phase.  All  boundary 
conditions  are  listed  in  Table  2. 

The  governing  equations  and  boundary  conditions  are  dis¬ 
cretized  by  the  finite  difference  method.  The  convergence  crite¬ 
ria  for  iteration  and  steady  state  are: 


ynew  T/old 

vi  vi 


v 


old 


<  1  x  10 


(12) 


V  is  an  arbitrary  variable. 


3.  Results  and  discussion 


The  parametric  study  for  a  two-phase,  PEMFC  cathode  tran¬ 
sient  model  is  presented.  The  parameters  used  in  this  study  are 
listed  in  Table  3. 

The  evolution  of  the  I-V  polarization  curves  are  shown  in 
Fig.  2.  When  the  cell  voltage  is  lower,  for  instance,  Vs  =  0.2  V, 
and  the  current  density  overshoots  with  time  due  to  the  rapid 
electrochemical  reactions.  But  when  the  cell  voltage  is  higher, 
for  instance,  Vs  =  0.6  V,  there  is  almost  no  overshoot. 

Fig.  3  shows  the  saturation  level  distribution  evolution  with 
Vs  =  0.2, 0.4,  and  0.6  V,  respectively,  since  there  is  no  water  in  the 


Table  3 

Parameters  used  in  simulation 


Gas  diffusion  layer  properties 

Porosity 

0.3 

Thickness 

0.025  cm 

Permeability 

10-9  cm2 

—(dpc/ds) 

284.2  Dyne  cm-2 

Catalyst  layer  properties 

Porosity 

0.06 

Thickness 

0.0016  cm 

Permeability 

3  x  10“ 11  cm2 

—(dpjds) 

568.4  Dyne  cm-2 

Catalyst  loading  (mpt) 

0.4 

Specific  surface  area  of  Pt  (apt) 

1000  cm2  (mg  Pt)- 

Volumetric  fraction  of  Nahon  in  catalyst  pellet  (e^) 

0.393 

Radius  of  catalyst  pellet  (rAgg) 

10-5  cm 

Thickness  of  Nafion 

10-6  cm 

Exchange  current  density  60, ref  0  °C) 

10-6  A  cm-2 

Membrane  properties 

Thickness 

0.005  cm 

Porosity 

0.35 

Fixed  charge  site  concentration 

1.2  x  10-3  mol  cm 

Operation  conditions 

Temperature 

60  °C 

Pressure 

1  atm 

Mole  fraction  of  oxygen  in  the  air  let 

0.206 

Humidity  in  air  inlet 

10% 

Humidity  at  anode 

100% 

gas  diffusion  layer  and  the  catalyst  layer  in  the  beginning.  After 
start-up,  liquid  water  is  generated  by  electrochemical  reactions 
in  the  catalyst  layer  and  then  diffuses  to  the  gas  diffusion  layer. 
The  liquid  water  saturation  level  accumulates  until  the  steady 
state  is  reached.  The  generated  liquid  water  will  begin  to  affect 
the  system  after  1  s.  Thus,  the  liquid  water  saturation  level  in  the 
catalyst  layer  is  higher  than  in  the  gas  diffusion  layer. 

When  the  cell  voltage  is  higher,  for  instance  Vs  =  0.6  V,  the 
electrochemical  reaction  is  moderate  and  less  liquid  water  is 
generated.  Thus,  the  liquid  water  saturation  level  is  lower  than 
when  a  higher  cell  voltage  is  applied.  When  the  cell  voltage  is 
lower,  more  liquid  water  is  generated  and  begins  to  occupy  the 
pores.  The  catalyst  pellets  surface  have  less  oxygen,  causing  a 
decreased  cell  performance  and  a  mass  transport  limitation. 


1  0.8 


1 .0 

Current  density.  [A  cm'2]  12  1.4 


Time,  t[sec] 


Fig.  2.  The  evolution  profile  of  polarization  curves. 
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Fig.  4.  The  evolution  profile  of  the  ionic  potential  £cl  =  0.1:  (a)  Vs  =  0.2  V,  (b) 
y -0.4  V,  (c)K=  0.6  V. 


Fig.  3.  The  evolution  profile  of  saturation  level  (base  case):  (a)  Vs  =0.2  V,  (b) 
y -0.4  V,  (c)K=  0.6  V. 
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Fig.  5.  The  evolution  profiles  of  current  density  with  various  gas  diffusion  layer  porosity:  (a)  Vs  =  0.2  V,  (b)  Vs  =  0.4  V,  (c)  Vs  =  0.6  V. 


Fig.  4  shows  the  evolution  of  the  ionic  potential  distribution 
with  Vs  =  0.2,  0.4,  and  0.6  V,  respectively.  In  the  membrane,  the 
ionic  potential  distribution  shape  is  a  straight  line  due  to  the 
absence  of  electrochemical  reactions.  In  the  catalyst  layer,  ionic 
potential  distribution  is  nonlinear  because  of  the  electrochemical 
reactions. 

When  the  cell  voltage  is  higher,  i.e.,  when  the  current  density 
is  smaller,  the  oxygen  reduction  reaction  and  the  ionic  potential 
loss  are  both  smaller  (Fig.  4c).  When  the  cell  voltage  is  lower, 
i.e.,  the  current  density  is  higher,  the  oxygen  reduction  reaction 
and  the  ionic  potential  loss  are  both  larger  (Fig.  4a).  The  present 
study  shows  several  phenomena  for  the  first  time. 

(1)  Firstly,  the  variation  of  the  ionic  potential  loss  reaches  a 
critical  value,  decreasing  to  a  steady  state,  and  is  not  mono¬ 
tonic. 

(2)  The  ionic  potential  reaches  80%  of  its  critical  value  in 
0.005  s. 

(3)  The  difference  between  the  critical  value  and  steady  state 
of  any  given  ionic  potential  is  dependent  on  cell  voltage 
applied.  When  cell  voltage  is  0.2  V,  the  difference  is  about 
20%. 


For  a  different  GDL  porosity,  the  cell  current  density  curves, 
with  Vs  =  0.2, 0.4, 0.6  V,  respectively,  are  shown  in  Fig.  5.  Within 
10-2  s,  the  current  density  rapidly  rises.  Between  10-2  and 
10“ 1  s,  the  current  density  remains  constant.  But  after  1  s,  the 


Fig.  6.  Effect  of  GDL  porosity  porosity  under  various  cell  voltages  in  the  steady 
state. 
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Time,  t  [sec] 
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Fig.  7.  The  evolution  profiles  of  current  density  with  various  catalyst  layer  porosity:  (a)  Vs  =  0.2  V,  (b)  Vs  =  0.4  V,  (c)  Vs  =  0.6  V. 


current  density  changes  again.  More  liquid  water  accumulates 
after  1  s,  as  shown  in  Fig.  2.  When  £gdl  =  0.2  and  £gdl  =  0.3,  the 
current  density  drops,  caused  by  liquid  water  hindrance  on  the 
oxygen  transport  due  to  reduced  GDL  porosity.  When  £gdl  =  0.4 
and  £gdl  =  0.5,  the  current  density  decreases  slightly.  This  is 
because  the  GDL  porosity  is  larger,  so  that  liquid  water  hindrance 
on  the  oxygen  transport  is  negligible.  When  Vs  =  0.6  V,  there  is 
no  liquid  water  effect  and  the  current  density  remains  constant. 
When  £gdl  =  0.5,  a  maximum  current  density  is  obtained. 

On  reaching  steady  state,  the  GDL  porosity  and  the  cell 
current  density  effects  are  shown  in  Fig.  6.  Under  the  same  cell 
voltage,  between  £gdl  =  0.2  and  £gdl  =  0.5,  gas  diffusion  layer 
porosity  increases  with  current  density.  When  the  cell  voltage 
is  lower  (Fig.  2),  the  liquid  water  effect  is  obvious.  When 
Vs  =  0.2  V,  the  current  density  difference  between  £gdl  =  0.2  and 
£gdl  =  0.5  is  more  than  40%,  and  when  Vs  =  0.6  V,  the  current 
density  difference  between  £gdl  =  0.2  and  £gdl  =  0.5  is  small. 

For  a  different  CL  porosity,  the  cell  current  density  curves, 
with  Vs  =  0.2,  0.4,  and  0.6  V,  respectively,  are  shown  in  Fig.  7. 
When  the  CL  porosity  is  higher,  the  catalyst  loading,  electro¬ 
chemical  reactions  and  the  current  density  all  decrease.  When 
£cl  >0.15,  the  current  density  is  lower.  But  a  too  small  CL  poros¬ 


ity  causes  the  liquid  water  effect,  the  same  as  for  lower  GDL 
porosity.  Furthermore,  when  £cl  =  0.04,  water  flooding  occurs 
even  at  10~2  s  and  when  Vs  =  0.6  V,  there  is  no  liquid  water  effect. 
The  current  density  increases  with  CL  porosity. 


Fig.  8.  Effect  of  CL  porosity  under  various  cell  voltages  in  the  steady  state. 
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On  achieving  steady  state,  the  CL  porosity  and  the  cell  current 
density  effect  are  shown  in  Fig.  8  from  which  it  is  clear  that,  when 
Vs  =  0.2,  0.4  V,  the  optimum  current  density  appears  between 
£cl  =  0.06  and  £cl  =  0. 10. 

4.  Conclusion 

A  parametric  two-phase,  one-dimensional  transient  model  of 
a  proton  exchange  membrane  cathode  is  the  basis  for  this  study. 
The  model  is  developed  to  investigate  the  transient  transport  of 
gaseous  species,  protons  and  liquid  water. 

Firstly,  the  variation  of  the  ionic  potential  loss  reaches  a  crit¬ 
ical  value,  decreasing  to  a  steady  state,  and  is  not  monotonic. 
This  phenomenon  is  particularly  obvious  at  lower  cell  voltages. 

The  current  density  rapidly  rises  within  10~2  s  and  remains 
constant  between  10-2  and  10-1  s.  But  after  1  s,  the  current  den¬ 
sity  is  affected  by  the  cell  voltage,  the  catalyst  layer  porosity,  and 
the  gas  diffusion  layer  porosity.  When  the  cell  voltage  is  higher, 
the  electrochemical  reactions  are  moderate,  with  no  liquid  water 
effects.  But  at  lower  cell  voltages  there  are  rapid  electrochemical 
reactions  and  a  liquid  water  effect.  When  Vs  <  0.4  V,  the  liquid 
water  effect  should  be  taken  into  consideration. 

For  the  gas  diffusion  layer  porosity,  the  current  density  drops 
when  £gdl  <  0-4.  This  is  caused  by  liquid  water  hindrance  on  the 
oxygen  transport  process  because  the  GDL  porosity  is  smaller. 
When  £gdl  =  0.5,  the  liquid  water  effect  is  small. 

For  catalyst  layer  porosity,  when  £cl<0.1  the  liquid  water 
effect  is  obvious,  but  when  £cl>0.1,  the  liquid  water  effect 
not  very  apparent.  When  there  is  not  enough  catalyst  to  partic¬ 
ipate  in  the  electrochemical  reactions,  the  current  density  goes 
down.  From  these  results,  an  optimum  value  appears  between 
£cl  =  0.06  and  £cl  =  0.1. 
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